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'^ ■ We study the effective interaction between a planar array of uniformly 

5 \ negatively charged, stiff rods parallel to a negatively charged planar substrate 

in the absence of salt in a continuous, isotropic dielectric medium. Using 
Brownian dynamics simulations, we examine the general effects of counterion 
valence, rod spacing, macroion charge densities, and the rod size on the at- 
tractive rod-surface interaction force. At room temperature divalent as well 
\f^ • as monovalent counterions mediate an interaction that can be repulsive or at- 

Tij- I tractive upon adjusting either the macroion charge densities or the rod radius. 

cn . Finally, we examine the effects of discretizing the surface charge as laterally 

^^ I mobile monovalent anions and of electrostatic images in the substrate. 

a^ 
a^ 

-(— > 

a 

o 
o 



X 



*Present address: NCSA, University of Ilhnois, 405 N. Mathews, Urbana, IL 61801; 
mashl@ncsa.uiuc.edu 
^ ngj@lanl.gov 



I. INTRODUCTION 

In the past few decades many investigations have provided evidence that effective at- 
tractions can occur between hke-charged macroions. AttractiorL between hke-charged sur- 
faces or rod-hke macroiatis^ has been observed numericallyailJ mainly through the use 
of Monte Carlo methods,ErLl but also through Brownian dynamics siniulationsly|-|Mpst of 
these studies have featured ions as mediating the attraction, but a few studiesQ'Ell have 
found an attraction by the use-ToL polyelectrolytes of the opposite charge. Attraction 
has also be^a^-seen analyticallyQ'EJn^ through-yieffiious methods, including integral-equation 
approaches,Qilic3-dpjisity-functional theories,!^" and lattice models.l^'E3 Finally, a variety 
of experiment g23i23E3 has indicated attractions as well. Foj^pcample, ^^^ay scattering stud- 
ies on colloidal particles made of poly(styrene sulfonate)E3c3 or silicaE3 have shown phase 
separation behavior, and when coupiad-iwith externaliv applied osmotic stress, has shown 
DNA aggregation by multivalent iong23E2 and ligands.E3 Effective diffusion constants derived 
from light scattering measurements on salt-free svstems of poly(styrene sulfonate) of varying 
molecular weights^ or of varying concentrationgi2l have suggested macroion bundling. Light 
scattering combined with electron microscopy has shown the reversible aggregation behav- 
ior of various filamentous biopolymers in the presence of multivalent metal ions or peptide 
homopolymers.cll Deuterium-NMR and ^^Na-NMR have been used to derive phase diagrams 
of lamellar liquid crystalline systems.c3 Finally, in calcium clays the attractive forces between 
the mica surfaces have been measured using the surface force apparatus,!^ and^the effects of 
osmolyte concentrations on the clays have been studied by X-ray scattering.ca Thus, there 
is substantial evidence of effective attractions between charged objects of the same sign. 

Theory has shown that ion-ion correlations (i.e., correlations in the fluctuations of the 
ion positions) are the source of the attractions in the numerically and analytically studied 
model systems. Mean-field descriptions, such as the Debye-Hiickel and Roission-Boltzmann 
equations,cj which are often used to calculate forces between macroions, EJlM^lways fail to 
yield this attraction because they do not account for interparticle correlations.!!^ On the other 
hand, integral equation approaches naturally possess the framework-for including such in- 
teractions. Fluctuation theories also can yield an attraction. Oosawa£3 has shown (to dipole 
order) that the coupling of longitudinal fluctuations in the counterionic atmospheres along 
two parallel, like-charged, rod-like macroions produces an electrostatic attraction. Whereas 
Oosawa has predicted this attraction to increase with temperature, we recently found bij 
the use of Brownian dynamics simulations that it increases as the temperature decreases.lSI 

Since the effective interaction between like-charged surfaces or two like-charged rods 
has been fairly well studied, we will address here the effective interaction between rod-like 
macroions and a like-charged surface. The geometry is motivated by experimental systems 
involving highly negatively charged DNA molecules and negatively chpiged membranes. In 
fact, experimental methods have been developed (for reviews, see Ref.E3EJ) to adsorb DNA 

se methods is the adhe^n of DNA to mica 
^ or cationic detergents.l^ The treatment of 



onto negatively charged substrates. Among-|t 
surfaces by the use of multivalent cations^ 
mica surfaces with a solution of magnesium acetate, calcium acetate, or magnesium chlo- 
ride has been seen to stabilKej-ajdsorbed DNA and subsequently givgiijeproducible imaging 
by atomic force microscopyE3^^ or scanning tunneling microscopy.E^TEll These treatments 
are believed to displace potassium ions in the mica with divalent cations so as to provide 
stronger binding sites. Similarly, divalent cations have been shown to be more successful 
than monovalent cations in binding DNA to quartz sand in DNA-degradation studies.E3 
Studies on gesje- delivery sy^ms have reported on the adsorption of DNA to anionic li- 
posomes withEj and without!^ the use of positively charged binding proteins, rather than 
simple cations. Finally, two-dimensional, fingerprint-like films of DNA strands Jbasie been 
found to adhere to cationic bilayers deposited onto untreated mica substrates.uC3 Thus, 
experiments reveal that adsorption of rod-like macroions onto like-charged substrates can 
be mediated by oppositely charged ions, proteins, or membranes. 

The model system we will investigate, as suggested the experiments above, is an array 



of rigid, negatively charged rods parallel to a negatively charged surface. The question on 
which we focus is, can a charged surface attract a like-charged rod? Although the DNA 
in the above experiments is admittedly flexible and the surfaces are not necessarily strictly 
planar, we investigate this simpler system in order to gain insight into the more complicated 
systems. For simplicity, we omit salt from our study in order to examine the fundamental 
nature of the interaction between the like-charged macroions. Our method is to fix the array 
of rods at certain distances from the surface, compute the time-averaged force on the rods 
perpendicular to the surface using Brownian dynamics simulations, and investigate the de- 
pendence of this force on a systematic variation of model parameters, such as array height, 
macroion charge densities, rod spacing, counterion valence, and rod size. We also examine 
the roles of electrostatic images and the discretization of the surface charge into mobile 
monovalent anions. In the next section we present the pairwise potential energy functions 
used in the simulations. These potentials model only the classical electrostatic interactions 
and neither distinguish between counterion species of the same valence nor consider solvent 
structure. We then summarize the Brownian dynamics algorithm by which the representa- 
tions of the counterions are moved by deterministic and thermal forces and follow with the 
results and a discussion. 



II. THEORY 



A. The Model 

The system we consider is an isotropic, continuous dielectric medium containing nega- 
tively charged rods and their counterions above a supported bilayer of anionic lipids. Our 
simulational unit cell model (Fig. 1) is a rectangular open half-space {z > 0) of dielectric 
constant ei (~ 80 for an aqueous medium) with dimensions L^ and Ly that contains a single 
line charge of uniform charge density A located a distance d from a fixed, fiat surface of 
a specified average charge density a. This surface density can be either a uniform den- 
sity distribution or a collection of discrete, mobile monovalent anions, representing, e.g., a 
fluid-phase mixture of charged and neutral lipids. The unit cell is replicated the x and y 
directions to produce an infinite array of infinitely long, parallel line charges with a spac- 
ing Lj, and with a repeat distance Ly in the y direction. There is no periodicity in the 
z direction. This replicated dielectric medium is bounded below {z < 0) by an isotropic, 
continuous half-space of dielectric constant €2 (~ 3) representing the dielectric layer formed 
by the hydrocarbon tails of the lipids. The dielectric discontinuity at z = produces elec- 
trostatic images in the substrate in-«rder to satisfy the boundary conditions on the dielectric 
displacement field at the interface.!^ A point particle of charge qi placed in the medium with 
ei at a distance Zi away from the planar interface will produce an electrostatic image charge 
of charge q'^ = gj(ei — e2)/(ei + £2) at the symmetric position —Zi. An electrostatic image 
rod is produced similarly. As a result, a uniform surface charge density (or the charge of 
individual mobile monovalent surface anions) increases by a factor of 2ei/(ei + £2). 

The interaction potentials consist of pairwise, long-ranged electrostatic (Coulombic) 
forces and short-ranged, non-electrostatic repulsions. We make no approximations regard- 
ing the Coulomb interations. For a uniformly charged surface, the ion-surface electrostatic 
interaction V^_s(z) is a function of the distance z the ion with charge qi is from the surface: 

vUz) - -^.. (1) 

where eo is the permittivity of vacuum. The rod-surface electrostatic contribution Vr_s(z) 
per unit cell for a uniform surface charge distribution is similarly 
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The ion- ion electrostatic potential energy l^_i(Ar) for an ion with charge gi at the 
point (x + Ax, y + A?/, z + Az) and an ion with charge g2 and its replicas located at 
(x + mLxiV + nLy,z), where m,n are integers, resulting from the two-dimensional repli- 
cation of the unit cell in the x and y directions, is given exactly by an absolutely convergent, 
double summation. The ion- ion potential energy l^_i(Ar) ig53 
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When the surface consists of discrete ions, the potentials Vi_s(-2) and K-s(-2) are identi- 
cally zero, and all electrostatic interactions with the surface are taken according to Eq. (^. 

The self-energy l^_/ that arises from an ion interacting with its own replicas is found by 
evaluating the expression 
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and is given in Ref.E3 The ion-rod potential energy V^_j.(Ar) is the logarithmic interaction 
between a point particle and a one-dimensional array of line charges. The analytic form 
of ^_r(Ar) is derived from the potential energy of an ion interacting logarithmically with 
a two- dimensional array of line charges arranged on a rectangular lattice (see Eq. (14) in 
Ref.O) by eliminating one of the dimensions. The result is given by 
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The corresponding ion-rod self-energy simplifies similarly from Ref.Ej to 
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The interaction between the real ions and electrostatic image ions as well as between the real 
ions (rods) and electrostatic image rods (ions) is obtained directly from the periodic poten- 
tials in Eqs. (0-§), respectively, up to a factor of 2ei/(ei -|- 62) ~ 1.93). Thus, interactions 
involving electrostatic images are treated without approximation. Finally, the ion-rod and 
ion-surface short-ranged, non-electrostatic repulsion energies were modeled as Ai_i./r^^ and 



Ai_s/r^°, respectively, to prevent the ions from collapsing electrostatically. The combination 
of the Coulomb interaction and short-ranged repulsion yields an optimal ion-rod distance, 
i.e., a rod radius tq, which was considered among the model parameters. Values for the 
coefficients Ai_j. and Ai_s and the resulting optimal ion-rod and ion-surface distances are 
given below. The ion-rod short-ranged interaction is taken according to the minimum image 
convention. 

The values of the parameters in our reference system are taken to be reasonably repre- 
sentative of thoscrlpund in the literature. Under physiological conditions the DNA found in 
most organisms isCi in the so-called B-form, a double-helical structure of approximate radius 
of 12 A, bearing two negatively charged phosphate groupsner base pair (i.e., every 3.4 A 
on average). Taken as a whole, calcium clays are reported" to have areas per unit charge 
ranging from 51-135 A^, whereas negatively charged lipids knownllHl to form (mono/bilayers) 

have a median surface area per molecule of about 51 A^. On the basis of this latter value 
alone, we choose our reference system to consist of rods with an equivalent linear charge 
density Arcf = — e/1.7A, lateral spacing L^ = 30 A, and a surface with an average charge 
density (Trcf = — e/5lA^. As we were unable to find data on the spacing of DNA in fully 
two-dimensional films adsorbed directly onto negatively charged substrate, we do not know 
whether a lateral spacing of 30 A is reasonab^ However, DNA molecules adsorbed onto 
cationic-lipid bilayers supported by mica haveEJ an interhelical spacing that depends non- 
monotonically on the concentration of divalent magnesium cations in solution and lies in the 
range 35-45 A. Furthermore, DNA molecules in sodium chloride solution are knownOo to 
pack into hexagonal arrays with lattice spacings ranging between 25-40 A. Dielectric con- 
stants in our system are taken to be ei = 80 and e2 = 3 or 80. In the former case, when 
electrostatic images are included, an ion or rod always repels its own electrostatic image. 
All simulations were carried out at 300 K. 

Values of A-r of 7.5 x 10^ 3.27 x 10^, and 1.8 x 10^ kcal A"/mol yielded rod radii 
ro of 5.6, 2.8, and 1.4 A, respectively, in the reference macroion system. Unless otherwise 
specified, the second Ai_r value was chosen as the reference rod value regardless of the actual 
macroion charge densities. For all simulations Ai_s was fixed at 6.63 x 10^ kcal A^°/mol, 
yielding an optimal ion-surface distance of 2.4 A for the reference system. The reference 
value of rod radius tq = 2.8 A is not equal to the real size of DNA. The reason for this is 
that the rod size and the correct ion-rod interaction energy cannot be imposed simultane- 
ously within our model. We believe it would be more correct to account for the electrostatic 
"binding" energy between a counterion and a DNA ion (on the orderE^ of a few ksT) rather 
than the actual DNA size. 



B. Simulation Method 

The motion of a collection of particles undergoing Brownian dynamics-|Corresponds to 
the long-time behavior of the Langevin equation in the overdamped limit.E3 The Langevin 
equation expresses a balance of deterministic and thermal, random forces on each particle 
at any given moment in time t. For each particle i with mass rrii, position rj(t), velocity 
Vj(t), and isotropic friction coefficient C,i, the Langevin equation reads 

m.^ = -^rmMt) + f.({r(t)}) + f;(t) (7) 

where f(({r(t)}) is the total deterministic force acting on particle i due to the positions 
{r(t)} of all the particles, and f*{t) is the random force acting on particle i. This random 
force is characterized by the fiuctuation-dissipation theorem as 



(C(0) ■ f;(t)) = 6m,kBn,5{t)5,, (8) 

for particles i and j, with the time average of the random force over any particle vanishing: 

(f;w) = o. (9) 

At long times the friction coefficent C,i is related to the diffusion constant D as Di = 
k-QT /mi^i. The inertial term on the left-hand side of Eq. (|^ can be set to zero at long 
times. One discretization of the resulting equation yields an iterative scheme for the time 
evolution of the system of particles:E3i3 

r,(t + At) = r,(t) + ^f,({r(t)}) + r*(At) (10) 

where r*(At) is a random Gaussian displacement with a variance of 2k-QT At / irii^i in each 
dimension. Given the total potential energy surface for a particle i, Vi_tot({r(t)}), the deter- 
ministic force fi({r(t)}) is given by fi({r(t)}) = -Vyi,tot({r(t)}). 



III. RESULTS AND DISCUSSION 

Each simulation required the specification of a set of system parameters {A, a, tq, L^., L^}, 
a fixed position of the rods d, whether the surface charge density consisted of discrete, mo- 
bile monovalent anions, and whether electrostatic images were to be included. The value 
of Ly depended on the particular set of system parameters and fell in the range 17-272 A, 
thereby giving approximately 80 or more ions in the basic unit cell to help ensure that 
particle configurations were well sampled. The divalent or monovalent counterions were 
introduced as a random configuration until the system became overall charge neutral; no 
additional ions (i.e., co-ions or salt) were added. As the rods were fixed in space, the simu- 
lation involved only the motion of the ions due to the explicit pairwise interactions among 
the ions, the rods, and the single charged surface, as well as the random thermal forces. 
For the ion motion, the normalized time step was taken to be At < 0.01 for m-j^j = 1. 
During the simulations the vertical and lateral time-averaged forces per unit length on the 
rod and the vertical time-averaged force per unit Ly length on the surface were monitored. 
After the initial 5-10 x 10^ steps were discarded, these averages were accumulated over an 
additional 2-200 x 10^ time steps, depending on the number of simulated particles, system 
size, macroion charge densities, and counterion valences. Equilibrium was attained when 
the average force on the rod (-F^ rod) and surface (-F^ surf) remained constant and the average 
lateral force on the rod (Fz^md) ~ 0. We checked for finite-size effects in Ly and found that 
they do not affect our results significantly. 

Numerous simulations with systematic variations of the system parameters were carried 
out to investigate the general behavior of the system. Figure 2 shows the time-averaged, 
rod-surface force in the 2;-direction per unit length rod as a function of rod-surface spacing 
for several rod and uniform surface charge densities for monovalent and divalent counterions, 
with and without electrostatic images. For the case of monovalent counterions, the top set 
of two curves shows that the effective rod-surface interaction is always repulsive, whereas the 
lower set of curves indicates that there are rod-surface distances d for which the interaction 
is attractive. As can be seen, the interaction becomes more attractive as the charge densi- 
ties on the rod and surface increase from one to four times the reference charge densities. 
For the case of divalent counterions, the rod-surface interaction is attractive for half, one 
times, and double the reference charge densities. The maximal (attractive) force, indicated 
by the depth of the force- distance profile, increases as the charge densities increase. Diva- 
lent counterions are obviously more effective at inducing an attraction than are monovalent 
counterions. 



The effect of changing the interred spacing L^ on a divalent-counterion system at the ref- 
erence rod and uniform surface charge densities was examined. For a system containing no 
electrostatic images, the force- distance profiles for L^ = 15, 30, 60 A were indistinguishable, 
and the profiles for systems with electrostatic images were also indistinguishable. Halving 
(doubling) the rod radius tq also resulted in indistinguishable profiles whose maximal attrac- 
tive forces, i.e., minima, were deeper (shallower) and located at smaller (larger) d compared 
to the reference case. Thus, these combinations of L^ and tq place these systems essen- 
tially in the dilute-rod limit. This dilute-rod regime is a well-defined limit as the effective 
rod-surface interaction is independent of the initial distribution of counterions. In further 
decreasing the lateral rod spacing L^ (and/or increasing the rod radius) one would even- 
tually encounter a less permeable array of rods. The resulting effective rod-surface force is 
then expected to depend longer on time and may lead to a significantly different effective 
rod-surface force as compared to the dilute-rod regime, due to the initial partitioning of 
counterions between the region exterior to the rods {z > d) and the region between rods 
and surface (0 < z < d, in Fig. 1). 

Figure 2 can also be interpreted in terms of the equilibrium rod-surface separation dis- 
tance, d*, the point at which force-distance profiles cross zero. In general, for either coun- 
terion species, as the macroion charge densities increase, the equilibrium rod-surface dis- 
tance decreases. This behavior is most striking for monovalent counterions as increasing the 
macroion charge densities causes d* to move from infinity to a finite distance. Finally, the 
effect of electrostatic images for both monovalent- and divalent-counterion systems is quite 
small compared to that which can be achieved through changing only the macroion charge 
densities. 

Changing the macroion charge densities individually leads to a type of phase diagram 
for the rod-surface force. Figure 3 shows the regions of rod-surface attraction and repulsion 
for systems of monovalent counterions with uniform charge densities and no electrostatic 
images. In general, the larger macroion charge densities (on the rod or the surface) lead 
to an effective attraction, and the lower charge densities lead to repulsion. The parameter 
set A/Aref = o ja-cci = 2 led to only marginal attraction, and this ambiguity is indicated in 
the figure. The dashed line, as suggested by the data, indicates only a qualitative phase 
boundary between the attractive and repulsive regions, as the actual shape of the boundary 
is unknown and would require a more detailed exploration of the macroion density phase 
space. Doubling the rod radius led to repulsion in each of the nine cases shown in Fig. 3. 
For a given fixed rod-surface separation distance, an increase in the rod radius tends to 
obstruct the counterions from entering the region in between the rod and surface, thereby 
hindering the ions from contributing relevant attractive interactions. Thus, the presence of 
an attraction bewteen the rods and surface is conditional on not only the macroion charges 
but also the rod radius. 

We examined the effect of macroion charge densities in systems of divalent counterions 
and found, as in Fig. 1, that an attraction between the rods and surface is evident in sys- 
tems with macroion uniform charge densities down to A = O.SArcf and a = O.Sarcf with the 
reference rod size. In these cases, the equilibrium rod-surface distance lies within a quite 
small range (5-6 A), with the smaller distances corresponding to the higher charge densities. 
We have not checked whether an attraction exists at even lower charge densities on either 
macroion. 

Finally, the roles of a discrete, mobile surface charge distribution and electrostatic im- 
ages in systems of monovalent and divalent counterions were examined. Figure 4 shows a 
comparison of systems with the reference densities A = Arcf and a = aj-et for the cases of 1) a 
uniformly charged surface with no electrostatic images, 2) a uniformly charged surface with 
electrostatic images, 3) a surface composed of discrete, mobile monovalent anions with no 
electrostatic images, and 4) same as in 3) but with electrostatic images. The force- distance 
profiles for the monovalent-counterion systems all show a repulsive rod-surface force and are 
essentially indistinguishable. Although the divalent force- distance profiles give nearly the 
same equilibrium rod-surface distances and appear similar in shape, the depths of the pro- 



files vary significantly. The maximal attractive force on the rods is greater for systems with 
surfaces composed of discrete, mobile charges than for systems with a uniformly charged 
surface. It is reasonable to attribute this result of enhanced attractions to the increased 
coupling between the mobile surface ions and the rod and other ions. A surface composed 
of a mobile charge distribution provides the system with more configurational degrees of 
freedom, which is, in general, expected to lead to an overall lowering of the free energy 
of the system. Interestingly, the four force profiles for the divalent-counterion systems all 
cross near 7 A. We do not know whether this common point is a general phenomenon or an 
artifact due to our particular system geometry and parameters. 



IV. CONCLUSIONS 

In summary, we have presented a systematic study of the effective interactions between 
a single charged surface and a parallel array of same-sign charged rods in the presence of 
their counterions with no co-ions (salt). We indeed find an effective attraction between the 
rods and surface that is mediated by either divalent or monovalent counterions at room tem- 
perature in the dilute- rod limit. This attraction results from counterion correlations that 
overcome the inherent repulsion between the two macroions. Divalent counterions readily 
mediate this attraction at macroion charge densities that are comparable to those found in 
highly charged biopolymers and surfaces. Clearly, counterions of higher valence are more 
effective at inducing attraction between charged HHjiacroions of the sanijs-sign charge density, 
and this trend has been observed experiment allyE^'El and numerically.Q'liy Although the very 
high macroion charge densities needed by systems of only monovalent counterions to display 
an effective attraction are seemingly unphysical, there indeed appears to be polymers with 
rather high equivalent linear charge densities. One groupE3 has pointed out that the anionic 
polymers dextran sulfate and heparin have a 2.0-fold greater charge density than DNA. The 
simulations have further revealed that the higher the macroion charge densities, the smaller 
the rod-surface equilibrium distance. 

We have demonstrated for this rod-surface system that the volume excluded by the rod 
to the counterions is an important factor in determining macroion attraction or repulsion, 
particularly for systems of monovalent counterions. An increase in rod radius tends to lower 
the attraction between the macroions because the average number of countericms in between 
the macroions is reduced. These results are consistent with our recent workll^l on systems 
of two identical rods in the absence of salt. In the present paper we also have discussed 
that the decrease in the maximal attraction due to an increased rod size can be partially 
offset by an increase in the rod-surface distance, while leading to a new decreased maximal 
attractive force. 

The accounting for electrostatic images generally introduces quite small, repulsive addi- 
tions to the rod-surface interaction. Whereas we have limited the counterions from coming 
within about 2.4 A from the surface, further investigation revealed that the inclusion of 
electrostatic images while decreasing the optimal ion-surface distance to about 0.5 A can in- 
crease the rod-surface equilibrium distance by 1 A compared to the case without electrostatic 
images. Thus, the inclusion of electrostatic images in these systems seems to introduce only 
rather short-ranged, repulsive effects. This raises a question for real aqueous systems and 
their modeling: To what extent do the hydration layers around macroions conceal electro- 
static image effects? 

There are many variations of this system that one could explore. For example, in the 
present system the effective rod-surface force could be studied as a function of the ratio of 
monovalent and divalent counterions. As the strengths of the electrostatic interactions in- 
volving monovalent counterions would be less than those involving divalent counterions, the 
"dilution" of a divalent-counterion system by the addition of monovalent counterions would 
be expected to decrease the attractive force between the rod and surface. Another variation 



of the present system geometry that has immediate relevance to experiments iSj-a system 
consisting of oppositely charged rods and surface. In particular, Yang and Fangp^ imaged 
DNA strands adsorbed onto a supported cationic lipid bilayer using atomic force microscopy 
and found that the DNA spacing in the fingerprint-like images increased monotonically from 
about 40 to 60 A as the sodium chloride concentration increased from about 20 mM to IM. 
This latter subject will be addressed in a future article. 
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FIGURES 




FIG. 1. The simulation system as described in the text. Periodic boundary conditions ap- 
phed to a unit cell of lateral dimensions Lx and Ly produce an infinite array of infinitely long line 
charges with a spacing of Lx- Each line charge with uniform charge density A and radius tq is 
located in an isotropic, continuous medium of dielectric constant ei and positioned a distance d 
from a charged surface at z = with average density a. The replicated ei medium is bounded by 
a second isotropic, continuous medium of dielectric constant 62- 
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FIG. 2. Force-distance profiles for the effective, timeaveraged force on the rod (per unit 
length) perpendicular to the surface for a given fixed rod-surface distance d. Shown are a series 
of curves for a system of monovalent counterions with (dashed lines, open symbols) and without 
(solid lines, filled symbols) electrostatic images. The symbols correspond to the following uniform 
charge densities: (•, o) A = Arcf , o' = cTvcf-, {^■, A) A = 2Arcf , cr = 2(Trcf ; (■, D) A = 4Arcf , o" = ^a^d- 
The inset shows the analogous series for a system of divalent counterions with the following charge 

densities: (•, o) A = O.SArcf , C = O.SiTrcf ; (A, A) A = Arcf, f = 0"ref ; (B, D) A = 2Arcf , O" = 2(Tref . The 

sizes of the error bars do not exceed those of the figure symbols and are not shown. In all cases, 
L^ = 30 A, ro w 2.8 A, Aref = -e/1.7A, and a^^i = -e/blh?. 



14 



n 



n 






Attractive 



n 



1 - + 



Repulsive 



I . \ . \ . \ . \ ^ 

12 3 4 5 

O / O r 

' ref 



FIG. 3. Phase diagram for a system of monovalent counterions, indicating regions in which 
the rod-surface interaction is attractive (□) and repulsive (+) for several values of uniform rod 
linear charge densities A and uniform surface charge densities a, relative to the reference system 
with Aref and Uref- The point at A/Arcf = o'/cref = 2 corresponds to a system that shows only 
marginal attraction (see Fig. 2) and is marked to indicate that it is a boundary point. The dashed 
line indicates qualitatively the boundary between the regions. In all cases, L^ = 30 A, ro ~ 2.8 A, 
Aref = — e/1.7 A and Uref = — e/51 A^. 
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FIG. 4. Comparison of systems with or without electrostatic images, having surfaces with ei- 
ther uniform charge densities or discrete mobile monovalent anionic charge distributions for A = Aref 
and a = iTrcf- The upper curves correspond to systems of monovalent counterions; the lower curves, 
divalent. The following key applies to both sets of curves: (+) uniform with no electrostatic images, 
(O) uniform with electrostatic images, (A) discrete mobile surface with no electrostatic images, 
(□) discrete mobile surface with electrostatic images. The sizes of the error bars do not exceed 
those of the figure symbols and are not shown. In all cases, Lx = 30 A, rg ~ 2.8 A, Afcf = — e/1.7A 
and (Tref = — e/51 A^. 
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